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Abstract 

We investigate the properties of the ground state of strong couphng lattice QCD at finite density. 
Our starting point is the effective Hamiltonian for color singlet objects, which looks at lowest order 
as an antiferromagnet, and describes meson physics with a fixed baryon number background. We 
concentrate on uniform baryon number backgrounds (with the same baryon number on all sites), for 
which the ground state was extracted in an earlier work, and calculate the dispersion relations of the 
excitations. Two types of Goldstone boson emerge. The first, antiferromagnetic spin waves, obey a 
linear dispersion relation. The others, ferromagnetic magnons, have energies that are quadratic in 
their momentum. These energies emerge only when fluctuations around the large- A'^c ground state 
are taken into account, along the lines of "order from disorder" in frustrated magnetic systems. 
Unlike other spectrum calculations in order from disorder, we employ the Euclidean path integral. 
For comparison, we also present a Hamiltonian calculation using a generalized Holstein-Primakoff 
transformation. The latter can only be constructed for a subset of the cases we consider. 
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I. INTRODUCTION 



The study of quantum chromodynamics at high density is of fundamental importance. 
In recent years the field has attracted wide attention with renewed interest in color super- 
conductivity (CSC). For a review see [1]. 

CSC at high density is a prediction of weak-coupling analysis. It is imperative to confirm 
this prediction by non-pcrturbative methods. Bringoltz and Svctitsky [2] studied high- 
density quark matter with strong-coupling lattice QCD (see also [3]). This theory is de- 
scribed by an effective Hamiltonian for color singlet objects. At lowest order in the inverse 
coupling, the Hamiltonian looks like an antiferromagnet, and describes the dynamics of 
meson-like objects within a fixed background of baryon number. Higher orders in the in- 
verse coupling add the effective Hamiltonian with kinetic terms for the baryons that are 
very hard to treat [2]. The global symmetry group of the antiferromagnet depends on the 
formulation of the lattice fermions. For naive, nearest-neighbor fermions, studied here and 
in Ref. 2, the symmetry is U{ANf). The representation of U{ANf) carried by the fields on 
each site depends on the baryon number on that site. 

From the antiferromagnetic effective Hamiltonian one can derive a path integral for a 
nonlinear sigma model (NLSM). The path integral can be studied in the large- A^'c limit with 
semiclassical methods. 

Following Rcf. 2, we choose to work with baryon number backgrounds that have the same 
baryon number B on all sites with \B\ — 1,2,..., 2A^/ — 1. These backgrounds are close to 
lattice saturation (that happens for B — 2Nf), but involve relatively simple calculations; 
The sigma fields on all sites belong to the same manifold and represent the same number 
of independent fluctuations (see [2]). Any other distribution of baryon number can also 
be treated with the NLSM, but will involve sigma fields whose structure changes from site 
to site. This makes the calculations of the semi-classical ground state and excitations less 
straight forward than presented here and in Ref. 2. We leave the treatment of non- uniform 
baryon number backgrounds, together with the kinetic terms for the baryons to future study, 
as a step towards comparing our results to those of grand canonical approaches. 

For the baryon number background we study, the classical {Nc = oo) vacuum of the 
NLSM is hugely degenerate: If one aligns the a fields on the even sites, the field on each 
odd site can independently wander a submanifold of the original a manifold. It was shown 
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in [2] that fluctuations in 0{1/Nc) couple the a fields on the odd sites to each other and 
give an effective action that hfts the classical degeneracy. The ground state breaks U (4A^/) . 

Here we treat the problem in the context of "order from disorder" [4-9], a phenomenon 
known in the study of condensed matter systems. Its essence is that fluctuations of quantum, 
thermal or even quenched nature can lift a classical degeneracy. A famous example is the 
Kagome antiferromagnet [7] whose classical ground state energy is invariant under correlated 
rotations of local groups of spins, leading to a degeneracy exponential in the volume. This 
system has modes with zero energy for all momentum. These zero modes obtain nonzero 
energy due to quantum fluctuations. The same thing happens in our NLSM. Here we identify 
the zero modes, show how they get nonzero energy of 0{1/Nc) and calculate their dispersion 
relations. 

We end up with two kinds of Goldstone boson. The flrst are antiferromagnetic spin 
waves with a linear dispersion relation. The second kind are ferromagnetic magnons with a 
quadratic dispersion relation. This is consistent with the loss of Lorentz invariance due to 
flnite density. This set of excitations falls into the pattern described by Chadha and Nielsen 
[10] and by Leutwyler [11] in their studies of nonrelativistic fleld theories. In general, a 
nonrelativistic system that undergoes spontaneously symmetry breakdown can posses two 
types of Goldstone boson. The energy of type I bosons is an odd power of their momentum. 
Their effective fleld theory has a second time derivative which means that each excitation is 
described by one real fleld as in the case of a relativistic scalar fleld. Type II bosons have 
a dispersion relation that contains an even power of the momentum. They can appear only 
in nonrelativistic theories having properties similar to a Schrodinger fleld theory, namely 
theories whose action has only a first time derivative. As in the Schrodinger case, each 
excitation is described by a complex field, or two real fields. The number of Goldstone 
fields is of course equal to the number uq of broken generators. The counting of massless 
excitations is summarized by the Chadha-Nielsen counting rule, 

ni + 2nn>nG, (1.1) 

Where ni{nii) are the number of type I and II Goldstone bosons. A well known example is 
the SU{2) spin system with a coUinear ground state , where uq — 2. The antiferromagnet 
has two spin waves with a linear dispersion relation {ni — 2, nu — 0) while the ferromagnet 
has one magnon with a quadratic dispersion relation (n/ = 0, nn = 1). As mentioned above 
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in our case both nj and rtn are nonzero. 

Most discussions of "order from disorder" work in a Hamiltonian (but see [8, 9]). We 
employ instead the Euchdean path integral. In Appendix C we also present a Hamiltonian 

calculation using a generalized Holstein-Primakoff [12] transformation. The latter can only 
be constructed for a subset of the cases we consider. The Euclidean calculations are less 
restrictive in their applicability. 

II. THE NONLINEAR SIGMA MODEL 

In Ref. 2 the problem of lattice QCD at finite densities was transformed to a nonlinear 
sigma model. In this section we review this approach in order to make the current work 
self-contained. The starting point is strong coupling QCD with naive fermions. This theory 
is described by an antiferromagnetic (AF) effective Hamiltonian, which is invariant under a 
global U{N) symmetry with N — ANf. This quantum Hamiltonian is given by 

-^eff= ^X)(3a/3(n)(5/3a(n + /i). (2.1) 

Here a and jS are Dirac-flavor indices taking values from 1 to N , and Qapi'n) are the gener- 
ators of U (N) on site n. The Hilbert space on each site forms an irreducible representation 
of U (N) which corresponds to a rectangular Young tableau with Nc columns. The number 
m of rows in the Young tableau can vary from site to site and is determined by the local 
baryon number according to 

m„ = 5„ + N/2. (2.2) 

Therefore, a first step in this Hamiltonian approach is to fix the local baryon number. This 
restricts the diagonalization of (2.1) to a certain sector of the Hilbert space, in contrast 
to other approaches that fix a chemical potential [1, 3] and let the local baryon number 
fluctuate. This point makes the comparison between the two approaches nontrivial and out 
of the scope of this paper. It also reduces the possibility of sketching a phase diagram to 
identifying the ground state for a given configuration of as was done in Ref. 2. 

To write the Euclidean path integral we use a generalized spin-coherent-state representa- 
tion, which turns out to be especially useful in the large- A^c limit. The result is a nonlinear 
sigma model (NLSM), which we proceed to describe. 
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A. Elementary fields 



Following [2] we study only the case of uniform baryon density, Bn — B > 0. The 
(7 field at site n is an N x N hermitian, unitary matrix that represents the coset space 
U{N)/[U{m) X U{N — m)]. It can be written as a unitary rotation of the standard matrix 

A, 

a„ = Ur^AUl (2.3) 

with 

—iN-m 

A convenient parameterization [12] for the unitary matrix Un is^ 

\ 

(2.5) 



A 



(2.4) 



1 - (Pn(Pn I 

Here 0„ is an m x (A^ — m) complex matrix that furnishes coordinates for U{N)/[U{m) x 



U{N — m)]. Inserting into Eq. (2.3), we obtain 

4 -1 + 2<f>l<Pr, 



\ 



(2.6) 



-2V1 

Note that Nc does not enter the definition of these degrees of freedom. 



B. The action and the ground state 



The action of the NLSM is 



2 



dr 



n ^ nix 



(2.7) 



The overall Nc factor allows a systematic treatment in orders of 1/Ac- The classical ground 
state is found by minimizing the action, which gives field configurations that are r indepen- 
dent and minimize the interaction. For B 0, this leaves a huge manifold of degenerate 
field configurations. 



^ This can be related to the form used in [2] by setting (j) = asin VcJa/ VaTa. The main advantage of this 
(j) parameterization is that the kinetic part of the action is quadratic in [see (3.6)]. 
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To clarify this important point, consider the energy of one hnk, 

^Tr (71(72. 



(2.8) 



To minimize E we work in a basis where ai = A. The analysis in [2] then shows that 
(72 can wander freely in the manifold U{m)/[U{2m — N) x U{N — m)], a submanifold of 
U (N) /[U (m) xU{N — m)] . A ground state of the infinite lattice can be constructed replicat- 
ing (7i and (72 on the even and odd sites of the lattice. Thus all the a fields on the even sites 
"point" to A while on the odd sites, each a field wanders independently in the submanifold. 
This classical ground state has a huge degeneracy, exponential in the volume. We proceed 
to analyze its stability. 

In [2] it was shown that 0{1/Nc) fluctuations lift the degeneracy and choose a coUinear 
ground state, where all the a fields on the odd sites point in the same direction. Without 
loss of generahty the ground state can be chosen to be 



= A, 



Codd 



A,= 








^ 

lAf-m / 



VAV\ 



(2.9) 



with 



1 



V 



(2.10) 



1 
0-10 

This ground state breaks U{N) to U{2m~N)xU{N-m)xU{N-m), with 2(3m-A^)(iV-m) 
broken generators. This was the main result of Ref. 2. 



III. FLUCTUATIONS AROUND THE GROUND STATE 

In this section we write an effective action for fiuctuations around the ground state up to 
0{1/Nc). Working at 0(1) first, we identify the zero modes that correspond to the classical 
degeneracy. We then show how quantum fiuctuations give them nonzero energy. 



6 



A. Effective action and 0(1) dispersion relations 



The fields 0, suitably shifted, can be identified with the Goldstone bosons around the 
ground state (2.9). For the even sites, = indeed gives a — K. For the odd sites, the 
vacuum is at 

^ 

where the upper part of is a (2m — N) x (A^ — m) matrix. We therefore shift on the odd 
sites according to 



u' = vu, 

(7'((/)) = Va{(l))VK 

This gives cr'(O) = A2. We drop the primes henceforth. 
For later convenience we write d) as 



(3.1) 







(3.2) 



Here tt is an {N — m)— dimensional square matrix while x has (2m — N) rows and {N — m) 
columns. Both are complex. Thus 

eleven = -27rX^ 1 - 27171^ -271-^1 - 0t0 , (3.3) 

^ -2^1 - 0t0;^t -2yi^0t^7r^ _i + 20V ^ 

^ l-2xxt -2x^1 - 0V 2x7rt ^ 
0-Odd = -2^1 - 0t0x^ _i + 20t0 2^1 - (P^n^ . (3.4) 

^ 27rx^ 27r^l - 0t0 1 - 27r7r^ ^ 

Each submatrix in Eqs. (3.3)-(3.4) has dimensions as indicated in (2.4). 
After scaling (p (p/y/l^ the action takes the form. 



and 



S ^ Tt dTj2 [xlidrXn + T^ndrT^r. 
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J 



n/j, 



_ t t _ t t 

'^n'^nXn+p.Xn+fi n+ pL^ n+ [iXriXn 

1 / t + 

+T^n+ixXn+fLXn+(,'^n + 7r„+A7r„xIiXri + /i-C 



+0(l/7Vf2)_ 



The AF structure of the action demands the introduction of an fee lattice. We write the fields 
0„ = 0;^ with A = (even, odd) and N belonging to an fee lattice. We Fourier transform 
according to 



^Ns^(3f-J' 6-^=^/2 A = odd 



(3.6) 



Here k = {k, cu). In momentum space the action is given by S — S2+S3+S4 with (discarding 
Umklapp terms) 



-^3 



-iu; + 2dJ -2 Jd7fc W \ 

Tr ^ ttI^ TT-i 

T ^ ^ \ -2Jd-fk ioj + 2Jd j 

= 2 J^^Tr {xfxlK-kli^ - x1xl-,K^H-r> + h.c) , (3.8) 

— 2 J(i-^^p— ^Tr ^ (XfeXp^XgXg+fe-p ~ ~ '^fe^''''p^g^g-fc+p 

5 ^' /cpq 



-X?X;«U-p - KKVxU+p) 7fe-P (3-9) 



7fc is given by 



1 

^fe = J Z)cosA;^/2. 



(3.10) 
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At large Nc, S3 and 5*4 are small perturbations. The bare propagators can be read from 



1 



^2 + AJ^d'^El 
I 



\ 



l/iu 
l/iuj j 



' iuj + 2Jd 2Jdjk 



AS 



(3.11) 
(3.12) 



AB 



(The propagators are diagonal in the internal group indices.) Here — y''l — 7^. The 
poles of the propagators give the dispersion relations of the various bosons at 0(1). The 
result is 



+ {2JdEkf = 



and 



to 







for TT, 



for X- 



(3.13) 



(3.14) 



The TT excitations are 2{N — rnf AF spin waves. A.i B — Q they are the only excitations. 
It is easy to verify that the Euclidean dispersion relation is a;^ ~ — jfcp at low momentum. 
The X excitations are A{N — m){2m — N) zero modes ("soft" modes). Their energy is zero for 
all momentum, a sign of the local degeneracy of the classical ground state discussed above. 

We conclude this subsection by classifying the different excitations according to their 
U {2m — N) X U{N — m) x U{N — m) representations. We denote a representation by 



(ri,r2,r3) 



(91:92,93) 



(3.15) 



Here ri denotes the representation of SU{2m — N), and r2,3 denote the representations of the 
two SU{N — m) subgroups, qi are the charges of the excitations under the remaining U{1) 
factors of the unbroken subgroup. These are generated by the following diagonal matrices. 

and 



i-N, 



i2m-N 






V 



(3.16) 





iN-m 
— lAr_m J 

We give the different representations in Table I. The crucial point is that the AF spin 
waves and the zero modes reside is completely different representations. This means that 
the separation between these excitations in the 0{1) spectrum will survive at higher orders 
and that mixing can not occur. 



Field 


Representation 


Dimension 


vr 


/ \ (0,0, +2) 

( 1, - m, - mj 


(AT - mf 


Xeven 


/ n(0,+1,+1) 

[2m-N,l,N-mj 


{2m-N){N-m) 


Xodd 


/ \(0,+l,-l) 


{2m-N){N-m) 



TABLE I: U{2m — N) x U (N — m) x U{N — m) representations of the excitations. The conjugate 
fields belong to the conjugate representations. n(n) stands for the fundamental(its conjugate) 
representation of SU (n) and 1 stands for a singlet. 



B. Self energy calculation and 0(1/Nc) dispersion relations of the x fields 



In this section we calculate the self energy of the x fields to first order in l/Nc- We shall 
see that the poles in the soft modes propagators move away from zero energy in this order. 
The Feynman rules for this problem are presented in Appendix A. We treat 5*3 and 5*4 as 
perturbations to 5*2 and present the contributions to the self energy in Fig. 1. 



.-(E>- = 



....0... 



+ 



+ 



FIG. 1: Dyson equation for the self energy of the x fields. Solid lines correspond to the vr propa- 
gators and dashed lines to the x propagators. 



Since the x propagator is divergent at a; = for all k, one must insert its self energy 
self- consistently, making the replacement 



E is a matrix coupling odd and even degrees of freedom according^ to 



E = 



El Ti2 
E* E* 



Thus the x propagator is replaced by 

1 



" (-ia;-Ei,fc)(ia;-E*,)-|E2,fe|2 



ioj - El J. E2,fe 



(3.17) 



(3.18) 



(3.19) 



We assume that the x propagator has a similar structure as Eq. (3.11). 
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Here we assume that both matrices Hi^k and T,2,k are diagonal in group space and that 

ReSi,fc<0 (3.20) 

for the stabihty of the path integral. 

In Appendix B we derive self-consistent equations for Ei^^ and Ti2,k- Defining tanh^fe = 
— Re E2,(o,fc)/Re Si (o,fc), we write these equations in the form of a single integral equation, 



hiq, k) sinh 9q 



(7V-m)tanh^fe= ^7 ^ = ^. (3.21) 

f]{k) and /i 2(q, fc) are defined in Appendix B. The poles of the propagators may then be 
obtained from (3.19) [see Eq. (B14)]. The dispersion relations turn out to be 

2Jd 



2Jd{N - m) 



2Jd{N-m) 



ak\l 1 - tanh^ Ok (3.22) 



The solution of Eq. (3.21) for d = 3 can be obtained numerically by assuming that the 
function 9k depends only on |fc| in most of the Brillouin zone. We plot the solution in Fig. 2 
for N — m = 1,2, 3, 4, 5. Recall that the baryon density is given by 5 = m — A^/2. Thus for 
Nf < 3 we cover all baryon density short of saturation. 

The form of efe is shown in Fig. 3. It is easy to check that the bosons are massless, since 
77(0) = -^1,2(9, 0) = implies that ccfe and Pk vanish at fc = 0. This is of course a direct result 
of the Ward identities concerning the global U (N) symmetry. Moreover, near k — both 
a and /3 have a quadratic dependence on |fe|. This means that at low momenta efc ~ |A; 
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characteristic of ferromagnetic magnons. A possible exception is the case N — m = 1. There 
we see that at low momenta, tanh6'fe — > 1 and thus it is possible that efe ~ Ifcp"*"^ with p a 
positive integer^. 

Finally we recall that the relation between the number of x fields and the number of 
physical excitations depends on the dispersion relation. Because the soft modes obey €k ~ 
the 4(A'" — m)(2m — N) fields describe only half that many ferromagnetic magnons. 
The case oi N — m — 1 might be different depending on whether p is odd or even. 



Cfe cannot be nonanalytic at A; = since the integrals in Eq. (3.21) are regular there. 
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|k| 



FIG. 2: Solution of Eq. (3.21) for - m = 1, 2, 3, 4, 5. 




FIG. 3: Rescaled energy of the x bosons for for iV — m = 1, 2, 3, 4, 5. The energy for iV — m = 1 is 
divided by a factor of 7 in order to present it on the same scale with the other cases. 

IV. SUMMARY AND DISCUSSION 



We find two kinds of low energy excitations corresponding to type I and type II Goldstone 
bosons [10, 11]. The type I bosons are 2{2Nf — BY antiferromagnetic spin waves that 
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have a linear dispersion relation. Including quantum fluctuations of 0{1/Nc) we also have 
4:B{2Nf — B) type II bosons^. We call them ferromagnetic magnons since they derive their 
energy from the effective ferromagnetic interaction [2] between next-nearest-neighbor sites, 
produced at 0{1/Nc). As typical for magnons, their dispersion relation is quadratic is 
momentum^ . 

The mechanism that removed the classical degeneracy and gave energy to the type II 
Goldstone bosons is known in the condensed matter literature as order from disorder [4-9] . 
It takes place in systems that posses a classical degenerate ground state. For example we 
mention the double exchange model [6] and the Kagome antiferromagnet [7]. There, the 
classical ground state energy is invariant under a rotation of local groups of spins making 
it exponentially degenerate, like in our problem. It is worth recalling that the zero density 
system is also classically degenerate [2, 13]. There, one can realize the baryon distribution 
by putting baryon numbers B and —B on the even and odd sites. At 0{1), the ground 
state energy does not depend on B. This leads to a discretely degenerate ground state. At 
0{1/Nc), order from disorder occurs and quantum fluctuation remove this degeneracy to 
pick out a ground state with B = 0. 

The energy scale of the type II Goldstone bosons is smaller by a factor of compared 
to that of the type I bosons. This points to a possible hierarchy of phase transitions at finite 
temperature which can be described by a classical model that has the Hamiltonian 

H = JiTvY^ [cr„a„+^] - J2 Tr ^ [crnan+2ii] ■ (4.1) 

np, nfi 

For J2 = 0, the ground state will be highly degenerate as above. A small positive value 
of J2 removes this degeneracy and picks out a ferromagnetic alignment of the next-nearest- 
neighbor a fields. Thus for T <^ Ji^2 the ground state is invariant only under 

U{N -m) xU{N -m) xU{2m- N). (4.2) 

This symmetry pattern will persists until some finite temperature Tj^^ ~ J2. For T > T™ 
the ferromagnetic magnons melt the magnetization and restore some of the symmetry. This 
situation will persists until a second temperature T^^^ where the symmetry will be 

^ Recall that the local baryon number B ranges from to 27V/ with unit increments. 

^ A possible exception is for N—m = 1, where the energy may depend on a different power of the momentum. 
Correspondingly, the number of excitations will be different. 
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restored completely. As long as Ji » J2 (which corresponds to A^c 1 in our system), the 
hierarchy of phase transitions is well defined. 

Finally we mention that recent works [14] on effective field theories for dense QCD also 
predict the existence of type II Goldstone bosons. It is tempting to identify these with our 
ferromagnetic magnons. First, however, one must reduce the artificial U{ANf) symmetry 
of the action to the physical chiral symmetry. We defer this important issue to further 
publications. 

APPENDIX A: FEYNMAN RULES 

The following are the Feynman rules extracted from Eqs. (3.7)-(3.10). A = (e, o) stands 
for even and odd respectively. Latin indices take values in the range [1, (A^ — m)^] while 
Greek indices take values in the range [1, (A^ — m)(2m — A^)]. The propagators are given in 
Fig. 4. The vertices extracted from the cubic and quartic interactions [Eqs. (3.8)-(3.10)] are 
given in Figs. 5- 6. 

(a) (Glm)^^ - ^ti ^ 4f iP) (Gl^^,,)^^ -xi. — 

FIG. 4: Propagators. 



= -2Jd^ Jfj^SapSjfdgi [dAo^Be^Co - (c ^ o)] 7fc. 

FIG. 5: Vertices corresponding to the cubic interaction. 
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A.ia 



^7/ 



XSg 




= "^Jdj^-J^SehSfiSjkSig [SAeSce^BoSDo + (c -(-^ o)] 7p_q, 




= ~'^d,jf^;^dih5jkSie5fg [SDe^AeSBoSce + (c •<-> o)] Jp. 



TT 



/9 




X = ^Jcg jy ^ g SgiSjeSjSsShf [SAeSceSBoSoo + (e ^ o)] 7p_ 



Xefj 



Xia 



TT, 





TT 



TT. 



tc 

e/ 



x; 



= -Jd j^l, g SasSjkSleSif [SAe^Be^Co^De + (e ^ o)] 7p. 



FIG. 6: Vertices corresponding to the quartic interaction. 
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APPENDIX B: THE SELF CONSISTENT EQUATIONS 



In this appendix we derive the integral equations represented by Fig. 1 and solve them 
approximately in order to obtain the dispersion relations of the x bosons to 0{1/Nc). In 
analogy with the propagators of the antiferromagnetic spin waves we assume the following 
structure for Ei^^ and I^2,k (we denote the frequency by fl), 



(Bl) 
(B2) 



where Ai_2, Bi^2 are real and Aj > A^ for stabihty. In order to perform frequency integrals 
we make a rational ansatz. For example, 



^ 2Jd P,in';k) 
2Jd P2{VL^\ k) 



(B3) 
(B4) 



N,Q2{n^;ky 

Where Pi_2(^^; k) and Qi,2(^^; k) are polynomials of the same order in a;^. Thus we have 



k)AB 



+ B^Y - Bl]+ Al - Al 



^m{l + Bi)-Ai A2 + iVtB2 ^ 
^ A2 - iVLB2 -iVL{l + Bl) - Ai J 



AB 



(B5) 



Since Bi^2 ~ 0{1/Nc) we drop them in the denominator. The self-consistent equations for A 
and B then decouple, and only the former are needed to find the poles of the x propagators. 
The equations® are 



-Pi 
Q 



El 



+ 



1 Pi{u:^-q)/Qi{u:^-q) 



{uj - Vtf + El N, u;2 + £;2(^2. 



7q - ll 



u:^ + E\uj^-q)\ [{uj - 9)^ + El_^ 
where the three terms come from the three diagrams in Fig. 1 and 



,(B6) 



P2(^'; fc) _ 2l 
Q2{n^; k) N, 



°° du; P2(oj^;q)/Q2(oj'^;q) 
J J -00 27r u;^ + E'^{u;^;q) 



X 



27fe7q - Iq-k 7 



7fc 



Iq-k 



(B7) 



® We have rescaled w 2Jdu) and Q. — > 2 JrfO. 
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which comes from the second and third diagrams in Fig. 1. Here 



= (B8) 

and 



p2 p2 IP 



It remains to evaluate the 0{1/Nc) contribution of the following integrals 

11 = P^^, (BIO) 
7^ = J- f % , (Bll) 

% . (B12) 

All these integrals are convergent at large cu and can be evaluated with residue calculus. 
Beginning with !{, 

The poles of the integrands are given by 

u;'Q'{u;') + ^P\u;')^0. (B14) 

c 

To leading order in 1/Nc the roots are determined by either 

= -^,P\0)/Q'{0) (B15) 

c 

or by 

g2(a;2) = 0. (B16) 

Since the polynomials Pi and Qi do not depend on Nc and Qi{0) can be chosen to be 1, the 
solutions of Eq. (B16) are all 0{1). Moreover since Q'^ > 0, all the roots appear in complex 
conjugate pairs u;*). For definiteness we choose Ima;„ > and close the contours from 

above. The leading contribution to Eq. (BIO) comes from the poles given by Eq. (B14). The 
calculation of I2 and I3 is similar, but one has two more poles to consider, at a; = O ± iE. 
The result is 

n = l^^, (BIT) 

2 7P2(0) _ p2(0) 
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r — 
-'2 — 



1 E 



(B18) 
(B19) 



29? + E^' 

Thus we find that to 0(l/7Vc), all integrals depend only on the values of the polynomials at 
a; = 0. Finally, Eqs. (B6)-(B7) simplify to 



(N-m) , , 

Jbz \47r^ 

E,-H {ll - it) 



cosh 1 + 



27fc7g7g-fc - 7q - 7r 



Where we have defined 



(dq 
BZ \ 47r 



sinh 6n 



EU + 



7q-fc 



tanh^, = ^i^. 



Taking 17 = 0, and dividing Eq. (B21) by Eq. (B20), we have 



{N - m) tanh Ok 

Here we have further defined 

r){k) = 
h{q,k) = 

l2{q,k) = 



oik 



( ^ 

BZ \ 47r 



d r 2 2 
Wq-fc - 7fc 



1 + 



27q7fc7q-fc - 7q - 7^ 
^ 27q7fc - 7q-fc (7q + 7fe) 



^q-k 



Iq-k 



A solution of this equation must obey 

/ 1 \ d 



BZ V An 



Ii{q, k) cosh 9q — rj{k) > 0, 



(B20) 
(B21) 

(B22) 



(B23) 



(B24) 
(B25) 

(B26) 



(B27) 



so that the stability condition (3.20) will be satisfied. 

After solving the integral equation (3.21) one can extract the spectrum of the x bosons 
by calculating the poles of the propagator (3.19). These poles are given by equation (B14). 
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As seen there are two kinds of solutions. The 0{1/N^ poles of the form (B15) give rise to 
a low energy band [see Eq. (B9)] 

2Jd 



= -^^4WM^P¥M (B28) 

Using Eq. (B22) and Eq. (3.21) we have 

2Jd(N-m) I -T— 2Jd(N-m) , 

= -J^ - tanh^ 9k = ^— Uk- (B29) 

The 0(1) poles that solve Eq. (B16) represent uninteresting massive excitations. 



APPENDIX C: HAMILTONIAN APPROACH FOR N-m = l 



In this appendix we derive the self-consistent equations in a Hamiltonian formulation 
using a generalized Holstein-Primakoff transformation ioT N — m — 1. 
The Hamiltonian is [2] 



J 



= H Qaf3{n)Qpa{n + ft). 



The U{N) generators Qapin) obey 



For N — m = 1 there is a simple representation of these operators [12], 

/ ljv-i-200t 20Jl-^t^\ 



Qal3 — 



\ 2^1 - 0t00t _1 + 20V ) 
Here is a (A'" — l)-component field that obeys 



(Cl) 



(C2) 



(C3) 



(C4) 



This representation is the generalized Holstein-Primakoff transformation. (We are not aware 
of such a representation for N — m> \. ) Next we write 







(C5) 
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The (A/" — 2) -component field x will represent the zero modes and the fields tt will constitute 
the antiferromagnetic spin waves. 

We write the generators on the even sites by inserting Eq. (C5) into Eq. (C3). As in 
Section III A we replace Q by 

g'(0) = VQ{<t^)V^ (C6) 

on the odd sites. V is given in Eq. (2.10). It is easy to see that the algebra (C2) on the odd 
sites is obeyed by Q' as well. 

Now, we scale (f) (j)/\fNl , substitute Q and Q' into Eq. (CI), and expand up to 
0(l/iVc). We are left with 



(C7) 



where 



H2 



n(i \ ij 

i 

+^ {T^nT^nT^i.T^n+fi + h.c.) - 7r„7rj:,7r„+/i7rJ,+^ + (n ^ Tl + /i)| 



(C8) 
(C9) 



(CIO) 



1. The ground state at 0(1) 



At lowest order in the Hamiltonian is given by H2, which does not depend on Xn- 

Moving to momentum space and diagonalizing H2 with a Bogoliubov transformation, we 
have 

N,H2 = 2Jd [iVe -{N- 1)] ^ ^Jl-ll (ala^ + b%) . (Cll) 
7q is given by Eq. (3.10) and 



^ cosh (fiq Sinh \ ( Ciq 

^ sinh (pq cosh cpq j y b_q 



(C12) 
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Here tanh 2(/7q = 7^ and the fields Oq and bq obey the usual commutation relations of ladder 
operators. The ground state |0) is thus given by 

aq|0) = bq\0) = 0. (C13) 

To 0(1) the excitations of |0) are AF spin waves with a linear dispersion relation. The ground 

state |0) has a local degeneracy, corresponding to arbitrary numbers of x bosons (creation 
of a X boson costs zero energy at this order). This is exactly the result of Section III A. 



2. Self energy and 0{1/N^ effective Hamiltonian 

In this section we calculate an effective Hamiltonian that splits the spectrum of the 
degenerate subspace discussed above. We use Rayleigh-Schrodinger perturbation theory 
[15]. H4 contributes at first order and i^s at second order. Both give a contribution of 
0{1/Nc). We must diagonalize the following Hamiltonian within the degenerate subspace. 

^eff = PH^P + PhJ^^H^P. (C14) 

P is a projection operator onto the degenerate sector and D is the energy denominator. Once 
^eflf is calculated, we follow [4] and use the Wick theorem to decouple its anharmonic terms 
in all possible ways. This step includes substitution of various bilinears by their vacuum 
expectation values. We calculate the vevs of the tt bihnears using (C12) and (C13), 

(TTTrt); = (TTTrt);^ - = |- sinh^ (C15) 

2 

and write an ansatz for the vevs of the % bilinears. 



= («±a) = ]^ E sinh ipq cosh ^qe^^'^^ (C16) 



{XiX% = {XiX%±(. = 1^ E ^1 A, (C17) 

(x£xi±A,^) = (x;.X;±Aj) ^-W^ ^2,S,je^"'^. (C18) 

Here Ai 2 are assumed to be real. Using Eqs. (C15)-(C18) we decouple H4 into the following 
form 

/ te 
Xki 



\X-ki 



(C19) 
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The calculation of of the second term in (C14) is very cumbersome. We will not pursue 
it here and just present the first steps. We take the projection operator to be 

l-P = Y.\ka){ka\ + \h){h\, (C20) 
k 

where, for example \ka) = flfclO) ® Here |x) is |fca)'s component in the x Fock space. 
In general (C20) is not correct, since it takes into account excitations of one spin wave only. 
Here it suffices since connects |0) only with excitations of that sort. Next we replace the 
energy denominator with 

D ^2JdN^^l--fl + 0(l). (C21) 

This assumes that the energy of the x bosons is smaller than the energy of the spin waves 
by a factor of N^.. After a lot of algebra we get 



k.i 



\ 




^ki 

\ X-ki J 



(C22) 



with 



/ (^] (Ai<, + l)/i(q,fc)-77(fc), (C23) 
Jbz \47ry 

h,2iq,k) and ri{k) are as given in Eqs. (B25)-(B26). To get self-consistent equations we 
calculate the vevs (C17)-(C18) using (C22). The resulting equations are 

+ , ^ (C25) 

^Jl - {J2k/Jlky 

A,, = , •^^'-/•^^'^ ■ (C26) 

Vi - {J2k/Jiky 

Defining Aifc + 1 = cosh^^, A2fe = sinh^fe and dividing Eq. (C26) with Eq. (C25) we get 
exactly Eq. (3.21) for - m = 1. 

This Hamiltonian calculation sheds light on the meaning of the function 9k in terms of 
vevs of the X operators. A stable sohition that obeys (3.20) has finite vevs. Mathematically, 
the integral equation (3.21) has also a trivial solution with 6'fc = that leads to Ai 2 = 0. 
However, this solution does not obey (3.20) and is thus unphysical, giving Ji < and making 
ifeflf not positive definite. 

Finally, we recall that the absence of a Holstein-Primakoff transformation foretold a 
Hamiltonian approach for A" — m > 1. The path integral approach is more general 
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